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ABSTRACT 

The large majority of extragalactic star cluster studies done to date have essentially 
used two or three-passband aperture photometry, combined with theoretical stellar 
population synthesis models, to obtain age, mass and extinction estimates, and some- 
• times also metallicities. The accuracy to which this can be done depends on the choice 

I of (broad-band) passband combination and, crucially, also on the actual wavelengths 

^SJ , and the wavelength range covered by the observations. Understanding the inherent 

^SJ ' systematic uncertainties (the main aim of this paper) is of the utmost importance for 

I a well-balanced interpretation of the properties of extragalactic star cluster systems. 

CO . We simultaneously obtain ages, metallicities and extinction values for ~ 300 clusters 

' in the nearby starburst galaxy NGC 3310, based on archival Hubble Space Telescope 

fH , observations from the ultraviolet (UV) to the near- infrared (NIR). We show that, for 

Q^' ages 6 < log(age/yr) < 9, and if one can only obtain partial coverage of the spectral 

energy distribution (SED), an optical passband combination of at least four filters 
including both blue and red passbands results in the most representative age distri- 
bution, as compared to the better constrained ages obtained from the full UV-NIR 
I SED coverage. We find that while blue-selected passband combinations lead to age 

distributions that are slightly biased towards younger ages due to the well-known 
age-metallicity degeneracy, red-dominated passband combinations should be avoided. 
k> I NGC 3310 underwent a (possibly extended) global burst of cluster formation 3 x 10^ 

■ yr ago. This coincides closely with the last tidal interaction or merger with a 

' low-metallicity galaxy that likely induced the formation of the large fraction of 

clusters with (significantly) subsolar metallicities. The logarithmic slope of the V- 
band cluster luminosity function, for clusters in the range 17.7 < F606W < 20.2 
mag, is apeoew — —1-8 ± 0.4. The observed cluster system has a median mass of 
(log(m/MQ)) ~ 5.25 ± 0.1, obtained from scaling the appropriate model SEDs for 
known masses to the observed cluster SEDs. 

Key words: Hii regions - galaxies: evolution - galaxies: individual: NGC 3310 - 
galaxies: interactions - galaxies: starburst - galaxies: star clusters 
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1 THE CLUSTER LUMINOSITY FUNCTION 
1.1 Introduction 

Cluster luminosity functions (CLFs) and colour distribu- 
tions are the most important diagnostics in the study of 
globular and compact star cluster populations in nearby 

* E-mail: grijs@ast.cam.ac.uk 



galaxies. For the old globular cluster (GC) systems in, e.g., 
the Galaxy, M31, M87, and old elliptical galaxies, the CLF 
shape is well-established: it is roughly Gaussian, with the 
peak or turnover magnitude at My ~ —7.4 mag and a Gaus- 
sian FWHM of ~ 3 mag (Whitmore et al. 1995, Harris 1996, 
2001, Ashman & Zepf 1998, Harris, Harris & McLaughlin 
1998). The well-studied young star cluster (YSC) popula- 
tion in the LMC and the Galactic open cluster system, on 
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the other hand, display a power-law CLF (Elson & Fall 1985, 
Harris & Pudritz 1994, Elmegreen & Efremov 1997). 

Hubble Space Telescope (HST) observations are con- 
tinuing to provide an ever increasing number of CLFs for 
compact YSC systems in galaxies beyond the Local Group 
and the Virgo cluster (e.g., Whitmore & Schweizer 1995, 
Schweizer et al. 1996, Miller et al. 1997, Zepf et al. 1999, de 
Grijs, O'Connell & Gallagher 2001, Whitmore et al. 2002). 
Although a large number of studies have attempted to detect 
a turn-over in young or intermediate-age CLFs, the shapes 
of such young and intermediate-age CLFs have thus far been 
consistent with power laws down to the observational com- 
pleteness thresholds (but see Miller et al. 1997, de Grijs et 
al. 2001, 2003a,b). 

1.2 Evolutionary Effects? 

The striking difference between the power-law distributions 

for young star clusters and the Gaussian distribution of the 
old Galactic GCs has recently attracted renewed theoreti- 
cal attention. The currently most popular (but, admittedly, 
speculative) GC formation models suggest that the distribu- 
tion of the initial cluster masses is closely approximated by 
a power law (e.g., Harris & Pudritz 1994, McLaughlin & Pu- 
dritz 1996, Elmegreen & Efremov 1997, Gnedin & Ostriker 
1997, Fall & Zhang 2001). 

Which processes will affect the CLFs such that they 
transform from a power-law shape to a Gaussian distribu- 
tion? It is generally assumed that the processes responsible 
for the depletion of a star cluster population over time-scales 
of a Hubble time include the preferential depletion of low- 
mass clusters both by evaporation due to two-body relax- 
ation and by tidal interactions with the gravitational field of 
their host galaxy (e.g.. Fall & Rees 1977, 1985, Elmegreen 
& Efremov 1997, MuraU & Weinberg 1997a,b,c, Ostriker & 
Gnedin 1997, Harris et al. 1998, Fall & Zhang 2001), and 
the preferential disruption of high-mass clusters by dynami- 
cal friction (Vesperini 2000, 2001). From the currently most 
popular GC evolution models it follows that any initial mass 
(or luminosity) distribution will shortly be transformed into 
a peaked distribution, although it should be noted that these 
models apply only to Milky Way-type conditions in which 
the GC system is characterised by significant radially depen- 
dent radial anisotropy. Vesperini (2000, 2001) has included 
the internal gravitational interactions between cluster stars 
in his models and concludes that these need considerable 
fine tuning to transform a power law initial cluster mass 
function (ICMF) into a Gaussian distribution, whereas a 
Gaussian ICMF conserves its shape rather independently of 
the choice of parameters: destruction of low-mass clusters by 
evaporation and the tidal field is balanced by the destruction 
of high-mass clusters through dynamical friction. 



1.3 Interacting galaxies 

AH of these models are valid only for Milky Way-type gravi- 
tational potentials, however. Galaxy-galaxy interactions will 
obviously have a major effect on the resulting gravitational 
potential, in which the dynamical star cluster evolution is 
likely significantly different (see, e.g., Boutloukos & Lamers 
2003, de Grijs et al. 2003a). 



The mass spectrum of molecular clouds or molecular 
cloud cores, progenitors of young star clusters, has never 
been determined in interacting galaxies to similarly faint 
limits as in the Milky Way and the Magellanic Clouds. The 
significant external pressure in interacting gas-rich galaxies 
may be expected to have a major effect on the molecular 
cloud mass spectrum. 

Significant age spread effects in young cluster systems 
- in which cluster formation is still ongoing - affect the ob- 
served CLF (Meurer 1995, Fritze-v. Alvensleben 1998, 1999, 
de Grijs et al. 2001, 2003a,b), which might in fact make an 
intrinsically Gaussian CLF (or ICMF) appear as a power- 
law CLF (see, e.g.. Miller et al. 1997, Fritze v. Alvensleben 
1998). It is obviously very important to age date the indi- 
vidual clusters and to correct the observed CLF to a com- 
mon age, before interpreting the CLF (Fritze-v. Alvensleben 
1999, de Grijs et al. 2001, 2003a,b), in particular if the age 
spread within a cluster system is a significant fraction of the 
system's mean age. 

In de Grijs et al. (2003a,b), we provide the first observa- 
tional evidence for a clear turn-over in the intermediate-age 
CLF of the clusters formed in the burst of cluster forma- 
tion in M82's fossil starburst region "B", which very closely 
matches the universal CLFs of old GC systems. This pro- 
vides an important test of cluster disruption models (see de 
Grijs et al. 2003b). 

1.4 Star cluster metallicities 

Metallicities of YSCs produced in galaxy interactions, merg- 
ers or starbursts are an important discriminator against GCs 
formed in the early Universe. They are expected to corre- 
spond to the interstellar medium (ISM) abundances of the 
interacting/starburst galaxies, and are therefore most likely 
significantly higher than those of halo GCs in the Milky Way 
and other galaxies with old GC systems. ISM abundances 
span a considerable range, however, among different galaxy 
types from Sa through Sd, irregular, and dwarf galaxies, and 
may exhibit significant radial gradients (Oey & Kennicutt 

1993, Zaritsky, Kennicutt & Huchra 1994, Richer & McCall 
1995). Many of these galaxies show abundance gradients, 
which have recently been shown to sometimes extend con- 
siderably beyond their optical radii (Ferguson et al. 1998, 
van Zee et al. 1998). However, from those large radii, gas can 
be funneled efficiently into the inner regions during galaxy 
interactions and mergers (e.g., Hibbard & Mihos 1995). 

Hence, a considerable metallicity range may be ex- 
pected for YSCs produced in interactions of various types of 
galaxies and even among the YSCs formed within one global 
galaxy- wide starburst. 

During a strong burst, typically lasting a few x 10* yr in 
a massive gas-rich galaxy, a significant increase of the ISM 
abundance may occur (Pritze-v. Alvensleben &: Gerhardt 

1994, their Fig. 12b). Meanwhile, some fraction of the gas 
enriched by dying first-generation burst stars may well be 
shock-compressed to cool fast enough to be built into later 
generations of stars or clusters produced in the burst. The 
same effect may occur when multiple bursts occur in a series 
of close encounters between two galaxies before their final 
merger. Hence, in extended starburst episodes metallicity 
differences between YSCs formed early on and late in the 
burst phase may be expected. 
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Precise (relative) metallicity determinations for individ- 
ual YSCs are not only important to study the questions 
raised above, but also for the correct derivation of ages from 
broad-band colours. 

Dust extinction is often very important in YSC systems. 
In particular the youngest post-burst galaxies and galax;ies 
with ongoing starbursts often show strong and patchy dust 
structures and morphologies. For instance, the youngest 
clusters in the overlap region of the two galactic discs in 
the Antennae galaxies arc completely obscured in the op- 
tical and can only be detected in near or mid-infrared ob- 
servations (Mirabel et al. 1998, Mengel et al. 2001). Older 
merger remnants like NGC 7252 or NGC 3921 seem to have 
blown their inner regions clear of all the gas left over from 
intense star formation during the burst (e.g., Schweizer et 
al. 1996). Extinction estimates towards individual YSCs are 
therefore as important as individual metallicity estimates in 
order to obtain reliable ages and to be able to derive an 
age- normalised CLF or YSC mass function. 

1.5 Spectroscopy vs. multi-passband photometry- 
Individual YSC spectroscopy, feasible today with 8m-class 
telescopes for the nearest systems, is very time-consuming, 
since observations of large numbers of clusters are required 

to obtain statistically significant results. Multi-passband 
imaging is a very interesting and useful alternative, as we 
will show below, in particular if it includes coverage of near- 
infrared (NIR) and/or ultraviolet (UV) wavelengths. The 
large majority of extragalactic star cluster studies done to 
date have essentially used two or three-passband aperture 
photometry, combined with theoretical stellar population 
synthesis models to obtain age estimates. The accuracy to 
which this can be done obviously depends on the number of 
different (broad-band) filters available as well as, crucially, 
on the actual wavelengths and wavelength range covered by 
the observations, and on the PSF size compared to the clus- 
ter surface density profile (i.e., on how close to the observa- 
tions were taken to the confusion limit for these clusters). 

In this paper we assess the systematic uncertainties in 
age, extinction and metallicity determinations for YSC sys- 
tems inherent to the use of broad-band, integrated colours. 
We have developed an evolutionary synthesis optimisation 
technique that can be applied to photometric measurements 
in a given rmmber N{N > 4) of broad-band passbands. 
The optimisation routine then simultaneously determines 
the best combination of age, extinction and metallicity from 
a comparison with the most up-to-date Gottingen simple 
stellar population (SSP) models (Schulz et al. 2002), to 
which we have added the contributions of an exhaustive 
set of gaseous emission lines and gaseous continuum emis- 
sion (Anders, Fritze-v. Alvensleben & de Grijs 2002, Anders 
& Fritze-v. Alvensleben 2003). We also compare these re- 
sults with similar determinations based on the Starburst99 
SSP models (Leitherer et al. 1999), but assuming fixed, so- 
lar metallicity for our sample clusters. Although this is an 
often-used assumption, we will show that this introduces sig- 
nificant systematic efi^ects in the final age distribution, and 
therefore also in the mass distribution. 

We decided to focus our efforts on the nearby, well- 
studied starburst galaxy NGC 3310, known to harbour large 
numbers of young star clusters, for which multi-passband ob- 
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Table 1. Overview of the HST observations of NGC 3310 





Exposure tiin6 
(sec) 


i^eiitre 




v_/r\,iEji\ J. 
(°) 


FSOOW 


900,1000 


WF3 


8645 


-150.282 


F336W 


2 X 500 


PC 


6639 


-144.232 


P439W 


2x300 


PC 


6639 


-144.232 


F606W 


500 


PC 


5479 


145.553 


F814W 


160,180 


PC 


6639 


-144.232 




100,160 


WF3 


8646 


-160.282 


FllOW 


2X159.96 


NIC2 


7268 


80.619 


F160W 


2X191.96 


NIC2 


7268 


80.619 



Notes: — Location of the galactic centre; ^ — HST programme identifier; 

— Orientation of tire images (taken from the image header), measured 
North through East with respect to the V3 axis (i.e., the X=Y diagonal of 
the WF3 CCD -H180°). 

servations from the near-UV to the NIR are readily available 
from the HST Data Archive (Section 2; see also Elmegreen et 
al. 2002, hereafter E02). In Section 4 we first place the NGC 
3310 starburst in its physical context. We then discuss the 
derived age distribution of the cluster population, which we 
extend compared to previous work, in terms of the evolution 
of the CLF and the interaction stage of its parent galaxy in 
Sections 5 and 6. We summarise our main results and con- 
clusions in Section 7. Finally, we will apply our knowledge of 
the systematic uncertainties gained in this paper to a larger 
sample of nearby starburst and interacting galaxies drawn 
from the HST GO-8645 programme (Windhorst et al. 2002) 
in Papers II and III (de Grijs et al., in prep.). 



2 OBSERVATIONS AND DATA 
PREPARATION 

NGC 3310 is a representative member of the class of galaxies 
often showing signs of active starbursts and recently formed 
star clusters. The galaxy may have been a normal, quiescent 
Sbc-type galaxy before it started to produce large numbers 
of new stars, possibly due to the merger with a gas-rich low- 
metallicity companion (cf. Kobulnicky & Skillman 1995). 

As part of HST programme GO-8645, we obtained ob- 
servations of NGC 3310 through the F300W ("UV") and 
F814W filters (Windhorst et al. 2002), with the galaxy cen- 
tre located on chip 3 of the Wide Field Planetary Cam- 
era 2 {WFPC2) . Observations in additional passbands were 
obtained from the HST Data Archive. In order to obtain 
the largest common field of view (FoV) in the optical wave- 
length range, we restricted these archival data to be taken 
with the WFPC2. In addition, we obtained archival NIR 
images taken with the Near-Infrared Camera and Multi- 
Object Spectrometer's (NICMOS) camera 2, which provides 
the best compromise of spatial resolution (0'.'075 pixel" ^) 
and FoV (19'.'2xl9'.'2) when combined with WFPC2 obser- 
vations (see below). We have summarized the combined set 
of HST observations through broad-band filters used in this 
paper in Table 1. 

Pipeline image reduction and calibration of the WFPC2 
images were done with standard procedures provided as part 
of the iraf/stsdas^ package, using the updated and cor- 

t The Image Reduction and Analysis Facility (IRAF) is dis- 
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rected on-orbit flat flclds and related reference files most 
appropriate for the observations. 

Wc first registered the individual images obtained for a 
given programme and created combined images for each pro- 
grammc/passband combination, using the appropriate IMA- 
LiGN and IMCOMBINE routines in iraf. Next, we rotated all 
of these combined images to a single orientation using ob- 
jects in common among the observed fields, applying iraf's 
ROTATE and IMGEOM routines. We used the observations of 
programme GO-6639 as basis for these rotations. Finally, we 
adjusted the pixel size of both the PC and the NICMOS/2 
images to that of the final, rotated WF3 images. The final 
pixel size for all images is 0'.'0998. 

Because of the significantly smaller FoV of the NICMOS 
observations, we created two sets of final, registered images: 
one set containing the WFPC2 images only, and a second 
set consisting of all aligned images {WFPC2 and NICMOS). 
The resulting combined "WFPC2-only" FoV is 331 x 323 
pixels (33'.'03 x 32'.'24). Similarly, the common FoV for all 
instrument/passband combinations is 226x227 pixels (22'.'55 
X 22'.'65). 

2.1 Source selection 

We based our initial selection of sources on a customized ver- 
sion of the DAOFIND task in the daophot software package 
(Stetson 1987), as running under idl. 

We obtained initial source lists for all available pass- 
bands, using a detection limit of four times the r.m.s. noise 
in the (global) sky background, determined from the indi- 
vidual images. We did not force our detection routine to 
constrain the source roundness or sharpness, in order to be 
as inclusive as possible. 

For NGC 3310, at a distance D ~ 13 Mpc (Bottinelli 
et al. 1984, Krcgel & Sancisi 2001), most star clusters and 
star forming regions appear as point-like sources, and the 
detection of individual stars among these is unlikely given 
the associated high luminosities required. Thus, we can use 
the initial source lists as a starting point to obtain the final, 
verified list of candidate star clusters, as outlined below. 

Wc cross-correlated the source lists obtained in each 
pair of subsequent passbands, allowing for only a 1.5-pixel 
positional mismatch between the individual images. The 
F606W (g) F814W cross-correlated source lists were subse- 
quently adopted as the basis for our final source lists, to 
which we added the source detections resulting from the 
cross correlations of the other passband pairs, if not already 
included. 

To reject artefacts remaining in the final source lists 
and real sources that were badly situated for aperture pho- 
tometry, we visually examined all of the candidates on en- 
largements of both of the images from which they originated 
for contrast, definition, aperture centering, and background 
sampling. We rejected candidates that were too diffuse or 

tributcd by the National Optical Astronomy Observatories, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National 
Science Foundation. STSDAS, tlic Space Telescope Science Data 
Analysis System, contains tasks complementary to the existing 
IRAF tasks. We used Version 2.2 (August 2000) for the data re- 
duction performed in this paper. 



that might be the effects of background fiuctuations. A small 
number of sources contained multiple components or sub- 
structure (as expected for young stellar associations), and 
the apertures were adjusted to include all of these. We also 
visually inspected the 2-D Gaussian fits to those sources that 
were assigned UGauss > 2 pixels in the initial automated fit- 
ting pass. If needed, apertures for our final photometric pass 
(and thus the source magnitudes), centre coordinates, and 
fGauss's were adjusted. The visual verification is a justifi- 
able step in the reduction process, since - after all - we are 
mainly interested in obtaining good photometry for a repre- 
sentative, sample of (young) star clusters, for which we can 
reliably derive the systematic uncertainties in their age, ex- 
tinction and metallicity determinations. Hence, we sacrifice 
some sample completeness in order to achieve higher sample 
reliability. 

Thus, our primary cluster candidate samples consist of 

well-defined sources with aoauss > 1-20 pixels and relatively 
smooth backgrounds. The total rmmber of visually verified 
sources contained in our source lists is 382 and 289, respec- 
tively, for the "WFPC2-only" FoV and for the combined 
FoV. 

2.2 Photometric calibration 

The coordinates from the source lists obtained in the previ- 
ous section were used as the centres for daophot aperture 
photometry in all passbands. 

The correct choice of source and background aperture 
sizes is critical for the quality of the resulting photometry. 
Due to the complex structure of the star-forming regions, we 
concluded that we had to assign apertures for source flux and 
background level determination individually to each cluster 
candidate by visual inspection. The "standard" apertures 
for the majority of the sources were set at a 5-pixel radius 
for the source aperture and an annulus between 5 and 8 pix- 
els for the background determination, although in individual 
cases we had to deviate significantly from these values (e.g., 
because of extended source sizes, contamination by nearby 
objects, or gradients in the local sky background). Our pho- 
tometry includes most, if not all, of the light of each cluster 
candidate. 

The photometric calibration, i.e., the conversion of 
the instrumental aperture magnitudes thus obtained to the 
-HST-fiight system (STMAG), was done by simply using the 
appropriate zero-point offsets for each of the individual pass- 
bands, after correcting the instrumental magnitudes for geo- 
metric distortions, charge transfer (in)eflBciency effects, and 
the exposure times, following procedures identical to those 
applied in de Grijs et al. (2002). 

The full data tables containing the integrated 
photometry of all clusters are available on-line, at 
http : //www . ast . cam . ac . uk/STELLARPOPS/Starbursts/. 

Special care needs to be taken when calibrating the 
UV aperture magnitudes, in particular the F300W and 
F336W fluxes. These fllters suffer from signiflcant "red 
leaks" (Biretta et al. 2000, chapter 3). This causes a frac- 
tion of an object's flux longward of ^ 6OOOA to be detected 
in these filters. Close inspection of the filter transmission 
curves reveals that the response curve of the red-leak re- 
gion of the F300W filter resembles the transmission curve of 
the F814W filter, being most dominant in the 7000-9000A 
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range. The red leak of the F336W filter is roughly simi- 
lar to the red half of the F814W filter response function. 
This is fortunate, and implies that we could in principle use 
our F814W observations to correct the F300W and F336W 
fluxes for red-leak contamination. Biretta et al. (2000) show 
that the red leak in these two filters is in general < 5% of 
the total flux for stellar populations dominated by K3V or 
earlier-type stars, although it can be as much as 10-15% 
of the total flux for stellar populations dominated by MOV- 
M8V stars. Thus, for starburst galaxies dominated by young 
hot stellar populations the red leak contamination should be 
almost negligible, typically < 1%. 

For realistic spectral energy distributions (SEDs) for 
early-type galaxies (dominated by older stellar populations) , 
we have shown that the F300W red leak is typically 5-7% of 
the total F300W flux, and never exceeds the 10% level, not 
even in the reddest galactic bulges (Windhorst et al. 2002, 
Eskridgo ot al. 2003). Adopting Eskridgc ot al.'s (2003) most 
straightforward assumption, namely that the red leak cannot 
account for any more than all of the observed counts in any 
given area of a few pixels in the F300W and/or F336W 
images, we derive a maximum contribution of the red leak 
in NGC 3310 to both the F300W and the F336W filters of 
3% of the total count rate in the F8I4 W image. We applied 
these corrections to the F300W and F336W images before 
obtaining the aperture photometry. 

2.3 Completeness analysis 

We estimated the completeness of our source lists by using 
synthetic source fields consisting of PSFs. We created arti- 
ficial source fields for input magnitudes between 18.0 and 
28.0 mag, in steps of 0.5 mag, independently for each of the 
F300W, F606W, and F814W passbands. We then applied 
the same source detection routines used for our science im- 
ages to the fields containing the combined galaxy image and 
the artificial sources. The results of this exercise, based on 
the WFPC2 FoV, are shown in Fig. 1. These completeness 
curves were corrected for the ofi^ects of blending or superpo- 
sition of multiple randomly placed artificial PSFs as well as 
for the superposition of artificial PSFs on genuine objects. 
A detailed description of the procedures employed to obtain 
these completeness curves was given in de Grijs et al. (2002). 

We found that the effects of image crowding are small: 
only < 1.5% of the simulated objects were not retrieved by 
the DAOFIND routine due to crowding. However, the effects 
of the bright and irregular background and dust lanes are 
large, resulting in variable completeness fractions across the 
central galaxy images. As a general rule, however, the curves 
in Fig. 1 show that incompleteness becomes severe (i.e. the 
corriplctenoss drops below ^ 50%) for F300W > 22.5 mag, 
for F606W > 24 mag, and for F814W > 24 mag. 

Foreground stars are not a source of confusion (e.g., 
E02). The standard Milky Way star count models (e.g., Rat- 
natunga & Bahcall 1985) predict roughly 1-2 foreground 
stars for the equivalent standard filter of F606W < 24 in 
the combined WFPC2 FoV of our images. 

Background objects may pose a (small) problem, how- 
ever, in particular among the fainter sources, since we did 
not impose any roundness or sharpness constraints on our 
extended source detections, in order not to omit unrelaxed 
young clusters from our final sample. However, such objects 




Magnitude 



Figure 1. Completeness curves for the WFPC2 FoV of NGC 
3310. The different line styles refer to different passbands, as in- 
dicated in the figure. 

should bo easily identifiable once we have obtained aper- 
ture photometry for our complete source lists, as they are 
expected to have significantly different colours. More quanti- 
tatively, background galaxies at redshifts greater than about 
0.1 are expected to have extremely red (F300W-F814W) 
colours compared to their local counterparts. 



3 ANALYSIS OF THE STAR CLUSTER 
PHOTOMETRY 

We applied a three-dimensional minimisation to the 
SEDs of our star cluster candidates to obtain the most likely 
combination of ago t, motallicity Z and total (i.e.. Galactic 
foreground and internal) extinction E{B — V) (assuming a 
Calzetti et al. [2000] starburst galaxy-type extinction law) 
for each individual object. Galactic foreground extinction is 
taken from Schlegel et al. (1998) for each individual object 
(Section 5; see Anders et al. 2002, 2003). The least-squares 
minimisation of the clusters' SEDs was done with respect 
to the Gottingen SSP models (Schulz et al. 2002), to which 
we added the age and metallicity-dependent contributions 
of an exhaustive set of gaseous emission lines and gaseous 
continuum emission (Anders et al. 2002, Anders & Fritze-v. 
Alvensleben 2003). Their contributions to the broad-band 
fluxes are most important for significantly subsolar abun- 
dances (i.e., 0.2Zq; because the YSCs and the ionised gas are 
assumed to have the same low metallicity) during the first 
~ 3 X lO'^ yr of their evolution (Anders et al. 2002, Anders & 
Fritze-v. Alvensleben 2003), while for solax metallicity pop- 
ulations the relative importance is significantly reduced (to 
roughly half the relative contribution of 0.2Zq objects) and 
of much shorter duration {t < 1.2 x 10^ yr). 

The Schulz et al. (2002) Gottingen SSP models compose 
a model grid in age and metallicity covering ranges of 0.02 < 
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Z/Zq < 2.5 (as discrete metallicities of 0.02, 0.2, 0.4, 1.0 and 
2.5Zq) and 1.4 x 10* < age/yr < 1.6 x 10^° yr, in time steps 
of 1.4 X 10* yr, respectively. They are based on the most 
recent Padova isochrones and include the thermally pulsing 
AGB stellar evolutionary phase for stars with masses m in 
the range 2 < m/M© < 7. We also extended the age range of 
the Gottingen SSP models down to ages as young as 4 x 10® 
yr in time steps as short as 4 x l(f yr for the youngest ages. 

In order to obtain useful fit results for all of our three 
free parameters, i.e., age, metalhcity and extinction"^" , we 
need a minimum SED coverage of four passbands (or three 
independent colours) . To assess the effects of choosing a par- 
ticular passband combination on the final results, we started 
with applying our fitting technique to passband combina- 
tions consisting of four filters each, and subsequently added 
more passbands, as follows: 

(i) a blue-selected passband combination, consisting of 
photometry in the F300W ("UV"), F336W {"U"), F439W 
( "B") and F606W ( "V") filters; 

(ii) a subset of only optical passbands: F336W, F439W, 
F606W and F814W ("/"); 

(iii) a red-selected passband combination: the red optical 
F606W and F814W filters, combined with the NIR FllOW 
( "J") and F160W ( "H") passbands; 

(iv) a five-passband optical/near- infrared combination: 
F439W, F606W, F814W, FllOW and F160W; 

(v) the previous combination with the addition of the 
blue F336W filter; and 

(vi) our full set of seven passbands from across the entire 
wavelength range. 

Finally, in Section 5.2 we compare our results from the 
Gottingen SSP models to independently determined solu- 
tions for the same star cluster samples based on the Star- 
burst99 SSP models (Leitherer et al. 1999) for ages below 
1 Gyr, and on the Bruzual & Chariot (2000; BCOO) models 
for older ages. The fitting of the observed cluster SEDs to 
the Starburst99 and BCOO models was done using a three- 
dimensional maximum likelihood method, "3/2DEF", with 
the initial mass Mi, age and extinction E(_B — V) as free 
parameters (sec Bik et al. 2003), and assuming a fixed, solar 
metallicity, a Salpeter-type initial mass function (IMF), and 
the (Galactic) extinction law of Scuderi et al. (1996). We 
will discuss the uncertainties introduced by the assumption 
of solar metallicity below. For the clusters with upper limits 
in one or more filters we used a two-dimensional maximum 
likelihood fit, using the extinction probability distribution 
for E(S — V). This distribution was derived for the clusters 
with well-defined SEDs over the full wavelength range (see 
Bik et al. 2003). 



t Strictly speaking, the cluster mass is also a free parameter. Our 
model SEDs arc valid for SSPs with masses of 1.6 X 10^ M©; to 
obtain the actual cluster mass, we scale the model SED to match 
the observed cluster SED using a single scale factor. This scale 
factor is then converted into a cluster mass. 



4 SETTING THE SCENE: THE YOUNG 
MERGER REMNANT NGC 3310 

NGC 3310 is a local, very active starburst galaxy with high 
global star formation rate, as is evident from its bright opti- 
cal and infrared continuum emission (Telesco & Gatley 1984 
[hereafter TG84], Braine et al. 1993, Mulder, van Driel & 
Braine 1995 [hereafter MvDB95], Diaz et al. 2000 [hereafter 
DOO], E02), strong and extended X-ray, UV, and thermal 
radio emission (van der Kruit & de Bruyn 1976, Code & 
Welch 1982, Fabbiano, Feigelson & Zamorani 1982, Zezas, 
Georgantopoulos & Ward 1998, Conselice et al. 2000), in- 
tense UV and optical emission lines typical of OB associa- 
tions (Heckman & BaUck 1980, Kinney et al. 1993, Mulder 
et al. 1995 and references therein) and its large global Ha 
equivalent width (MvDB95, Mulder & van Driel 1996 [here- 
after MvD96], and references therein). 

4.1 A peculiar starburst galaxy 

A number of morphological peculiarities in its outer parts 
(e.g., Balick & Heckman 1981, Kregel & Sancisi 2001, [here- 
after BH81, KSOl], MvD96), combined with the disturbed 
kinematics of the Hi gas in the inner regions (e.g., van der 
Kruit 1976, MvDB95, KSOl), the mismatch between stel- 
lar and gas-dynamical geometry (e.g., BH81), and the large 
number of early-type stars required to explain the galaxy's 
Hq emission (van der Kruit 1976, BH81), suggest that NGC 
3310 was affected by a major gravitational disturbance. This 
led to high, possibly sustained, star formation rates in the 
past ~ 100 Myr (cf. BH81). Since attempts to identify a 
nearby companion galaxy as the cause for the disruption and 
the expected subsequent major starburst were unsuccessful 
(cf. van der Kruit 1976, van der Kruit & de Bruyn 1976), it 
was suggested that NGC 3310 accreted a gas-rich, but metal- 
poor companion galaxy, which subsequently fragmented as 
a result of the encounter (BH81, Schweizer & Seitzer 1988, 
MvDB95, MvD96, Smith et al. 1996, hereafter S96; modeled 
by Athanassoula 1992 and Finer, Stone & Teuben 1995), or 
perhaps we are currently seeing a newly-formed disc. This 
argument was based predominantly on the absence of any 
close companion galaxy and on the unusually low (subsolar) 
metallicity found in star forming regions surrounding the 
nucleus, although the nucleus itself appears to have solar 
metallicity (e.g., Heckman & Balick 1980, Puxley, Hawar- 
den & Mountain 1990, Pastoriza et al. 1993, hereafter P93). 
Additional support for this interpretation is provided by the 
far-ultraviolet and S-band R^^'^ luminosity profiles, typical 
of late-stage galaxy mergers (S96, KSOl), and the optical 
and Hi tails and ripples observed in the outer parts. 

The galaxy's most conspicuous visual peculiarity is 
the well-studied bright optical "bow and arrow" structure 
(nomenclature first used by Walker & Chincarini 1967). The 
"arrow" is a jet-like structure with possibly an Hi coun- 
terpart (although the latter is significantly more extended; 
MvDB95), which has been interpreted as the result of a nu- 
clear explosive event some 1.8 x lO'^ yr ago (Bertola & Sharp 
1984), or - in combination with the ripple pattern that in- 
cludes the "bow" - as the result of a recent merger with 
a smaller companion galaxy (or even a third, smaller ob- 
ject; KSOl). The latter explanation seems more likely (see 
MvDB95, KSOl). These authors argued that a "projectile" 
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object hitting the disc of NGC 3310 at a fairly obhque angle 
could have been drawn out into the observed configuration, 
which is supported by the anomalous velocity structure ob- 
served within these features (MvDB95); the optical "arrow" 
is then best understood as a scries of regions in which active 
star formation was triggered due to the collision. If this sce- 
nario is correct, the merger event must have occurred < 10 
Myr ago to prevent significamt redistribution of the longer 
Hi jet due to difi^erential rotation (cf. MvDB95), although 
the outer jet regions appear to curve away from the line of 
sight (KSOl, see also MvDB95). 

4.2 Star cluster formation 

The bar-driven star formation scenario suggested by Con- 
selice et al. (2000), combined with the recent infall of a com- 
panion gala^xy, is very attractive: it provides a natural expla- 
nation for the low metallicity observed in the star-forming 
knots near the centre, while it also explains why we observe 
concentrated star formation in star clusters or very luminous 
Hll regions in a tightly-wound ring-like structure surround- 
ing the centre (e.g., van der Kruit & de Bruyn 1976, TG84, 
P93, Meurer et al. 1995, S96, Consehce et al. 2000, E02), 
coinciding with the end of the nuclear bar (Conselice et al. 
2000, but see DOO), but not inside this ring. 

The most luminous single star-forming region in NGC 
3310 is the giant "Jumbo" Hii region (BH81, TG84, S96); 
it is roughly 10 x more luminous than 30 Doradus in Ha 
(BH81), but it is of very low metaUicity, Z ~ O.IZ3 (P93). 
The Jumbo Hll region contains several individual UV-bright 
star clusters (Meurer et al. 1995, E02). The most luminous 
star-forming structure is the ring of HlI regions, which pro- 
duces some 30% of the observed far-UV flux of NGC 3310 
(S96). We have indicated our cluster detections in Fig. 2, 
overlaid on our press release image of the galaxy (Wind- 
horst et al. 2001). 

Numerous other, less luminous star clusters and Hii re- 
gions, with sizes from 10 pc for the most luminous ones 
to < 1 pc for the unresolved clusters in the background 
disc (see Conselice et al. 2000), are found scattered in the 
galaxy's outer parts beyond Ji ~ 4 kpc (van der Kruit & de 
Bruyn 1976, BH81, TG84, MvDB95, S96) and as conden- 
sations in the "bow and arrow" structure (e.g., Bertola & 
Sharp 1984, S96, KSOl); these clusters may have been pro- 
duced by the accretion event, or might be remnants of the 
progenitor companion galaxy (cf. BH81, S96). 

4.3 Star formation time-scales 

All of the observational evidence points at very recent star 
formation in the star clusters and Hll regions, and a time 
since the interaction of < 10^ — 10^ yr (van der Kruit 1976, 
BH81, TG84, P93, S96, E02): the optical colours and the 
equivalent widths of the Ha-bright circumnuclear sources 
are best reproduced by a combination of a 2.5 Myr and an 
8 Myr-old population (cf. P93, DOO). This is consistent with 
(i) the detection of WR features in the spectra of a few of 
the Hii regions, including the Jumbo region, and of the NIR 
Ca H triplet (Terlevich et al. 1990, P93), both indicative 
of stellar populations with ages < 4 Myr, and (ii) the ab- 
sence of significant non-thermal radio emission (cf. TG84), 



which implies that there is not yet a significant population of 
supernova remnants. Terlevich et al. (1990) found evidence 
for two stellar populations (of 5 and 15 Myr old) as well, 
which can also be interpreted as an extended star formation 
episode. As pointed out by TG84 (see also MvDB95), the 
star formation in NGC 3310 can continue at its present rate 
only for a small fraction of the age of the galaxy, roughly 
4 X 10^ - 10* yr, depending on the IMF assumed (MvDB95). 



5 AGE DETERMINATIONS: SYSTEMATIC 
UNCERTAINTIES 

5.1 Dependence on the SED coverage 

We obtained 156, 145, 100, 85, 63 and 60 solutions from our 
least-squares minimisation technique (for our age, metallic- 
ity and extinction determinations), respectively, for the six 
passband combinations chosen for the analysis of the NGC 
3310 photometry (see Sect. 3) in the small FoV covered by 
the full passband combination. The relatively small fraction 
of acceptable solutions mainly reflects the rather low S/N 
ratios in a number of our archival images, resulting in ei- 
ther large photometric uncertainties or upper limits for the 
aperture photometry. Since the main aim of this paper is to 
analyse the dependence on the choice of passband combina- 
tion of the resulting ages, extinction values and metallicities, 
objects containing upper limits to their photometry in any 
of the passbands in a given passband combination were ex- 
cluded from the fitting procedure. We emphasize once again 
that we used exactly the same sample of star clusters, with 
well-determined photometric measurements in all passbands 
(as defined above), for this comparison. Acceptable solutions 
were required to have fit values for each photometric mea- 
surement in a given passband combination within the 2a ob- 
servational uncertainty associated with that particular filter; 
weights were allocated proportional to 1/x^- 

Fig. 3 shows the dependence of the resulting age distri- 
bution on the choice of passband combination. The shaded 
histograms represent the acceptable age solutions (with min- 
imum values obtained from the least-squares minimisa- 
tion) for the cluster sample in common among all six pass- 
band combinations; the open histograms are the additional 
solutions obtained for the respective passband combinations. 
The uncertainties in the number of clusters in a given age 
bin are predominantly Poissonian; the uncertainties in the 
actual ages are the subject of our current analysis. The 
electronic data tables containing our best solutions for all 
passband combinations are also available from the afore- 
mentioned WWW address. It is immediately clear from this 
figure that the resulting age distribution of an extragalactic 
star cluster system based on broad-band colours is a sensi- 
tive function of the passbands covered by the observations. 

To guide the eye, we have overplotted - on all panels 
- the best-fit Gaussian age distribution based on the full, 
7-passband SED sampling (Fig. 3f). This best-fit age distri- 
bution is characterized by a peak at log(age/yr) = 7.51±0.12 
and a Gaussian width of crcauss = 0.33, irrespective of the 
bin size adopted. 

From a comparison of the age distributions in Fig. 3, 
we determine that: 

• the peak of the age distribution is robustly reproduced 
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Figure 2. Cluster detections in NGC 3310, for wliicli wc obtained acceptable solutions for all of their ages, mctallicitics, and extinction 
values (see Sect. 5.1) overlaid on a grey-scale rendition of our F300W press release image of the galaxy (Windhorst ct al. 2001, 2002). 
The circles represent sources younger than log(age/yr) = 7.1, the squares are the older objects. The nomenclature of sources A-C, E, 
and L is from P93. 



in all cases, for both the shaded and the full cluster sam- 
ples with acceptable solutions, except for the red-selected 
passband combination (Fig. 3c); 

• photometric measurements in red-dominated passband 
combinations (Figs. 3c and d) result in significantly older 
(but highly uncertain) age solutions. This is due to the weak 
time dependence of the NIR magnitudes and due to un- 
avoidable ambiguities in the modelling of the thermally puls- 
ing AGB phase, thus resulting in badly constrained "best" 
fit ages. Red-dominated age solutions produce a significant 
wing of older clusters, compared to the ages derived from 
the full passband combination. 

• blue-selected passband combinations (such as in Figs. 
3a and b; open histograms) tend to result in age estimates 
that are slightly skewed towards younger ages, compared to 
passband combinations that also include redder passbands. 
This is simply due to the combination of observational selec- 



tion effects, in the sense that younger clusters have higher 
UV and [/-band fluxes and are therefore more easily de- 
tected at those blue wavelenghts, and to the age-metallicity 
degeneracy (see below) . Such secondary peaks at young ages 
should therefore be treated with extreme caution. 

• the optical passband combination (Fig. 3b), in particu- 
lar, displays a distinct secondary peak at very young ages, 
log(age/yr) ~ 6.7 (cf. open histogram). These clusters do 
not appear to be exceptional in their metallicity, extinction 
or mass properties compared to the overall distribution of 
these properties among the other clusters in our sample. 
However, with few exceptions, we could not obtain satisfac- 
tory solutions for these clusters as part of our full passband 
combination fits. As we will show below, this secondary peak 
is most likely the result of the age-metallicity degeneracy. 

In Fig. 4 we illustrate the sensitivity of the parameter 
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Figure 3. Relative age distributions of the NGC 3310 clusters based on six choices of broad-band passband combinations, with acceptable 
solutions. The shaded histograms represent the cluster sample with acceptable solutions in common among all six passband combinations 
(60 clusters); the open histograms are the additional solutions obtained for subsamples based on the passband combination displayed 
in each panel. The numbers in between brackets in each panel correspond to the total number of clusters in the histograms. We have 
overplotted the best-fit Gaussian age distribution obtained in panel f. 



determinations on the available passband combination by 
showing the SEDs and best fits for three typical, representa- 
tive clusters, located in P93 regions "C" (Fig. 4a; cluster No. 
18) and "B" (Figs. 4b and c; Nos. 116 and 122, respectively). 
The dotted lines represent the best-fitting models based on 
a blue-selected passband combination, the dash-dotted lines 
were derived from the red-selected combination, while the 
dashed and solid fits are based on the full passband combi- 
nation and on the subset of UBVI filters only, respectively. 
In Table 2 we list the best-fitting free parameters as a func- 
tion of passband combination used for all three clusters, to 
illustrate the effects of our choice of wavelength coverage. 

Note that all of our conclusions apply to the YSC popu- 
lation in NGC 3310, characterised by ages in the range from 
~ 10® to ~ 10® yr. The situation may be significantly dif- 
ferent for clusters of significantly greater age. Thus, from an 
observational point of view, for ages 6 < log(age/yr) < 9, 
we conclude that if one can only obtain partial coverage of a 
star cluster's SED, an optical passband combination includ- 
ing both blue and red optical passbands results in the best 
balanced and most representative age distribution, as com- 
pared to the better constrained ages obtained from the full 
UV-NIR SED coverage. While blue-selected passband com- 
binations lead to age distributions that are slightly biased 
towards younger ages, red-selected passband combinations 
- in particular if they are dominated by NIR filters - should 
clearly be avoided. 



Table 2. Dependence of our fit parameters on the available pass- 
band combination for three representative clusters 



Passbands log( Age ) log{ Mass ) Mctallicity E{S - V) 

[yr] [Mq] (Z) (mag) 



Cluster No. 18 (from Region "C") 



UV-!7SV 


7.205 


5.666 


0.022 


0.05 


UBVI 


6.876 


5.356 


0.008 


0.00 


VI JH 


6.660 


5.304 


0.013 


0.25 


W-UBVIJH 


7.204 


5.512 


0.002 


0.00 




Cluster No. 


116 (from Region 


"B") 




W-VBV 


7.362 


5.968 


0.011 


0.08 


UBVI 


7.448 


5.939 


0.008 


0.06 


VI JH 


8.939 


6.770 


0.002 


0.06 


W-UBVIJH 


7.478 


5.911 


0.012 


0.06 




Cluster No. 


122 (from Region 


"B") 




UV-UBV 


6.986 


5.474 


0.012 


0.13 


UBVI 


7.276 


5.577 


0.010 


0.10 


VI JH 


8.071 


6.079 


0.003 


0.12 


VV-UBVIJH 


7.349 


5.636 


0.012 


0.11 



5.2 The age— metallicity degeneracy 

Finally, in Fig. 5 we compare the best-fit age distributions 
of the clusters in the small FoV covered by all seven pass- 
bands, and the larger "WFPC2-only" FoV, covered by our 
five optical passbands. These best-fit age distributions were 
derived from either the full, seven/five passband SED cov- 
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Figure 4. Illustration of the sensitivity of the parameter deter- 
minations on the available passband combination for three rep- 
resentative clusters in regions "C" (panel a) and "B" (panels b 
and c). The dotted lines represent the best- fitting models based 
on a blue-selected passband combination, the dash-dotted lines 
were derived from the red-selected combination, while the dashed 
and solid fits are based on the full passband combination and on 
the subset of UBVI filters only, respectively. The HST filters are 
spaced according to their central wavelengths. 

erage if acceptable solutions could be obtained for all of our 
free parameters, or from the passband combination with the 
highest "confidence ranking", in a minimum sense: we 
defined a confidence ranking among the passband combi- 
nations, such that a higher ranking was given to solutions 
obtained from a passband combination (i) containing more 
passbands, (ii) covering a larger wavelength base line and 
(iii) resulting in a smaller value. 

For comparison, in Fig. 5a we have - once more - over- 
plotted the best-fit Gaussian age distribution obtained from 
the full SED sampling of the same (small) FoV, as in Fig. 3f. 
The vertical error bars represent the Poissonian uncertain- 
ties dominating the cluster numbers in each age bin. In Fig. 
5b we show both the age distribution of the clusters in the 
large FoV (solid histogram) and again the age distribution 
in the small FoV, overplotted as the dashed histogram (both 
normalised by the total number of clusters in each sample) , 
in order to illustrate the apparent ~ (2— 2.5)cr peak at young 
ages, log(age/yr) < 6.8. 

A naive interpretation of this feature would suggest that 
the outer field of NGC 3310, outside the small FoV, contains 
a subpopulation of significantly younger clusters. However, 
closer inspection of our data reveals that this young peak 
is entirely owing to unavoidable systematic uncertainties in- 



herent to our fitting technique and the shape of the SED. 
It turns out that the secondary, young peak is caused by 
the limited number of passbands available for the full clus- 
ter sample in the WFPC2-only FoV, due to the lack of NIR 
photometry for the clusters outside the small FoV. If we 
restrict our fits of the age, extinction and metallicity prop- 
erties of the cluster sample in the small FoV to the same 
WFPC2-only passbands, we find a similar young age peak 
in the inner FoV, of similar strength and significance. 

The clusters contained in this latter young-age peak are 
characterized by a smoother age distribution (ranging from 
~ 5 X 10® to ~ 2 X 10^ yr, with a few outliers as old as 
~ 10* yr) in the best-fit age distribution based on the full 
UV-NIR SED sampling. In addition, while these clusters are 
all characterized by a narrow spread in metallicity around 
~ O.5Z0 in the restricted fit using the WFPC2 passbands 
only, their metallicity estimates range from ~ 0.02Zq to 
~ 2.5Z0 if we take the fuU UV-NIR SED into account. 
Thus, we conclude that the secondary age peak at young 
ages in the WFPC2-only FoV is an artefact caused by the 
age-met allicity degeneracy at ages younger than about 10* 
yr. The age-extinction degeneracy is not as significant as 
this age-metallicity degeneracy, for these particular clusters. 
The slight, ~ 2cr excess of young clusters in the youngest 
age bin of Fig. 5a may be a remaining systematic effect, 
or may indeed be marginal evidence of more recent cluster 
formation. 

To illustrate this point, in the top panel of Fig. 6 we 
show the age dependence of the model SEDs for solar metal- 
licity. Next, in the bottom panel we explore the effects of 
varying the metallicity for the youngest, 8 Myr-old SED. 
For illustrative purposes we also show the effects of adding 
E(_B — V) = 0.5 mag extinction to the solar-metallicity SED 
(dashed line with solid bullets). The grey vertical lines are 
drawn at the central wavelengths of the F606W (or F555W) 
and F814W HST filters often used for age dating of individ- 
ual clusters. It is clear that for these young ages, the effects 
of varying the metallicity, extinction or age of the stellar 
population cannot be disentangled unless one has access to 
observations in a larger number of passbands, which may 
allow to more robustly constrain all of the free fitting pa- 
rameters. 

Finally, in Fig. 7 we compare our best fit age distri- 
bution (dashed Gaussian distribution) with the best fit age 
distribution using the Bik et al. (2003) approach. The latter 
is based on the Starburst99/BC00 SSP models, and a fixed 
metallicity of Z = Zq. While the differences between the 
age distributions resulting from fitting the Gottingen or the 
Starburst99/BC00 models are expected to be random due 
to the use of different stellar evolutionary tracks and dif- 
ferent template spectra for (super)giant stars (see Schulz et 
al. 2002), the assumption of fixed, solar metallicity causes 
the fitting technique to produce a significantly difi'erent age 
distribution compared to the distribution obtained by leav- 
ing the metallicity as a free parameter. We confirmed this 
result by fitting the individual cluster ages and extinction 
values using the Gottingen SSP models, while keeping the 
metallicity fixed at solar values. This is, again, a clear sig- 
nature of the age-metallicity degeneracy. We point out that 
this is, in fact, an expected effect, since unusually low (sub- 
solar) metallicities were found in NGC 3310's star forming 
regions (e.g., Heckman & Balick 1980, Puxley, Hawarden & 



© 2002 RAS, MNRAS 000, 000-000 



Star Clusters in Starburst Galaxies. I. 11 



3 



0.4 



0.3 



0.2 



k 0.1 



(a) Small FoV 



7 



0.4 



0.3 



0,2 



0.1 



(b) Large FoV 



log( Age / yr ) 



Figure 5. Comparison of the age distributions of the NGC 3310 star clusters in the small, inner FoV and the large, WFPC2-only FoV. 
The error bars represent Poissonian uncertainties in the numbers of clusters in each age bin. The dashed Gaussian distribution (panel 
a) is our best- fit age distribution (Fig. 3f); the dashed histogram overplotted on the large FoV age distribution (panel b) is the age 
distribution of the small FoV. 



Mountain 1990, P93, DOG; Z ~ (0.2 - 0.4)^0), so that the 
assumption of solar metallicity for these sources must have 
a significant systematic effect. 

5.3 Metallicity and extinction estimates 

Figure 8 shows our best estimates of the global distribu- 
tion of metallicity, and of total extinction (i.e.. Galactic 
foreground and internal extinction) towards the NGC 3310 
clusters. At first sight, the metallicity and extinction dis- 
tributions obtained from either of the passband combina- 
tions used for Fig. 8 appear very similar within the (sys- 
tematic) uncertainties. The effects of the age-metallicity 
degeneracy for the optical passbands (Fig. 8a) are clearly 
visible. Based on the UB VI passband combination only, we 
find that the metallicity distribution of the clusters in NGC 
3310 is strongly dominated by significantly subsolar metal- 
licities, while the addition of NIR passbands results in the 
detection of a significant number of clusters with metallic- 
ities Z ~ 2.5Zq. In either case, the fact that the result- 
ing metallicity determinations are strongly dominated by 
(significantly) subsolar metallicities is encouraging, in view 
of independent metallicity measurements in the literature. 
Both the low and higher-metallicity objects are distributed 
fairly smoothly across the face of the galaxy; there is some 
evidence that the most actively star forming regions, in par- 
ticular the Jumbo region and the northern spiral arm (in- 
cluding complex C), are predominantly composed of lower- 
abundance star clusters. 

Similarly, the addition of NIR passbands appears to re- 
sult in more solutions with lower extinction estimates (Figs. 
8e and f). This is likely due to the better constraints on 
the extinction from optical-NIR compared to optical-only 
baselines. However, whether this result is indeed signifi- 
cant also depends on the suitability of the starburst galaxy- 



type extinction law adopted for our fitting routine (Calzetti 
et al. 2000), of which the analysis is beyond the scope of 
the present paper (but see Calzetti, Kinney, & Storchi- 
Bergmann [1994] for a discussion). Nevertheless, these ef- 
fects are negligible, or marginally significant at worst, com- 
pared to the photometric uncertainties (~ 0.10 — 0.15 mag). 

6 NGC 3310 IN ITS PHYSICAL CONTEXT 

6.1 Overall properties of the NGC 3310 cluster 
population 

We can conclude from this exercise that NGC 3310 under- 
went a significant burst of cluster formation some 3 x lO'^ 
yr ago. The actual duration of the burst of cluster forma- 
tion may have been shorter because uncertainties in the age 
determinations may have broadened the peak. It appears, 
therefore, that the peak of cluster formation in NGC 3310 
coincides closely with the suspected galactic cannibalism or 
last tidal interaction (van der Kruit 1976, BH81, TG84, P93, 
S96, E02), while the possibly more recent, marginally sig- 
nificant cluster formation at t < l(f yr can be interpreted 
as cluster formation associated with the most recent (< 10 
Myr) enhanced star formation episode traced by the circum- 
nuclear Ha-bright sources (P93, DOO). 

The clusters older than ~ 10^ yr are smoothly dis- 
tributed throughout the galactic centre, roughly following 
the inner ring structure and the other concentrations of star 
clusters (see Fig. 2). However, the younger clusters, with 
ages log(age/yr) < 7.1 (i.e., the youngest two age bins in 
Fig. 3), are predominantly concentrated in the Jumbo re- 
gion ( "A" ; nomenclature from P93) and in the northern spi- 
ral arm, which also contains the bright star-forming regions 
C and E (see Fig. 2). This is consistent with the very young 
ages derived for these regions, t < 10 Myr (Terlevich et 
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Figure 8. Distributions of metallicity, and of total extinction towards the NGC 3310 clusters for the solutions from a number of passband 
combinations. The numbers in between parentheses correspond to the total number of clusters represented by each histogram. 



al. 1990, P93, DOO, E02), although wo emphasize that both 
regions also contain clusters spanning the entire age range 
observed for the NGC 3310 star cluster system. 

Figure 9 shows our best estimates of the clusters' 
masses. These were obtained by scaling our model SEDs 
(for masses of 1.6 x 10^ M©) to the observed SEDs for the 
appropriate combination of ago, metallicity and (total) ex- 
tinction. Our model SEDs arc based on a Salpeter-type IMF 
consisting of stellar masses m in the range 0.15 < tti/Mq < 
(50 — 70) (Schulz et al. 2002). The exact upper mass limit 
for our model SSPs depends on the metallicity and is deter- 
mined by the mass coverage of the Padova isochroncs. The 
exact value for the upper mass limit is unimportant for the 
determination of the total cluster masses. 

We realise that recent determinations of the stellar 
IMF deviate significantly from a Salpeter-type IMF at low 
masses, in the sense that the low-mass stellar IMF may well 
be significantly flatter than the Salpeter slope. The impli- 
cation of using a Salpeter-type IMF for our cluster mass 
determinations is therefore that we may have overestimated 
the individual cluster masses (although the relative mass dis- 
tribution of our entire cluster sample remains unaffected). 
Therefore, we used the more modern IMF parametrisation 
of Kroupa, Tout & Gilmore (1993, hereafter KTG) to de- 
termine the correction factor, C, between our masses and 
the masses obtained from the KTG IMF (both normalised 
at I.OM0). This IMF is characterised by slopes of a = —2.7 
for m > I.OM0, a = -2.2 for 0.5 < rn/Me < 1.0, and 
-1.85 <a< -0.70 for 0.08 < m/Me < 0.5. Depending on 
the adopted slope for the lowest mass range, we have there- 
fore overestimated our individual cluster masses by a factor 



of 1.5 < C < 2.4 for an IMF containing stellar masses in the 
range 0.15 < ui/Mq < 70. 

The mass distributions derived from the different pass- 
band combinations are fairly similar. We determined the 
defining parameters of the cluster mass distributions by 
fitting Gaussian distributions to them. The small differ- 
ences among the median mass [ (log (m/M©)) = 5.22,5.29, 
and 5.21 for Figs. 9a, b and c, respectively] and Gaus- 
sian widths of the distributions (similarly varying among 
CGauss = 0.32, 0.34, and 0.34, respectively) represent the 
minimum systematic uncertainties in these parameter deter- 
minations (minimum because they axe based on our highest- 
confidence passband combinations). In view of the uncer- 
tainties introduced by the badly known lower-mass slope of 
the adopted IMF, the median mass estimates of the NGC 
3310 cluster system are remarkably close to those of the 
Galactic GC system and of the intermediate-age cluster sys- 
tem in M82 B (de Grijs et al. 2003a,b). 



6.2 Comparison of individual cluster results with 
previous determinations 

Although both NGC 3310 as such and its star cluster system 
in particular have been studied extensively, determinations 
of individual cluster ages, masses, metallicities and extinc- 
tion values are rare. The most extensive cluster sample use- 
ful for a comparison to the results presented in this paper 
was published by E02. They estimate masses and ages of (a 
subset of) 11 large-scale, diffuse cluster complexes (see also 
DOO) and of 17 super star cluster (SSC) candidates based 
on HST photometry; the latter were also noted by Conselice 
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Figure 9. Mass distributions for the NGC 3310 clusters for a number of passband combinations. Overplotted are the best-fit Gaussian 
distributions (soUd lines). The numbers in between brackets correspond to the total number of clusters represented by each histogram. 



et al. (2000). In Table 3 we compare the E02 mass and age 
determinations for their sample of SSC candidates, based on 
the assumption of a mean metallicity oi Z = 0.008 = OAZq, 
with the corresponding parameters derived in this paper. 

We give our mass estimates for two assumptions of the 
lower mass limit of the Salpeter-type IMF; varying the up- 
per mass limit has very little effect. While we have used a 
lower mass limit of 0.15Mq throughout this paper (resulting 
in the mass estimates Mtot,i), E02's cluster mass estimates 
(Mtot.cf ) are based on the Starburst99 models, which assume 
a lower mass cut-off of IMq. Everything else being equal, 
our resulting mass estimates are therefore expected to be a 
factor of ^ 13 higher compared to theirs. Our total mass 
estimates assuming a lower mass cut-off of IMq are referred 
to as Aftot,2. 

A first comparison of our individually derived metallic- 
ity estimates shows that E02's assumption of a mean clus- 
ter metallicity of OAZq (based on spectroscopic metallic- 
ity determinations by P93) is reasonable, except for a few 
clusters with significantly supersolar abundances. Secondly, 
our low extinction values are also in line with our expecta- 
tions for the circumnuclear ring clusters; although the ex- 
tinction in the starburst ring varies in the range 1 < Ay < 4 
mag (Grotheus & Schmidt-Kaler 1991), is has been shown 
that circumnuclear ring clusters are often either almost fully 
obscured or virtually dust free (e.g., Maoz et al. 2001). In 
addition, Meurer et al. (1995) estimate the average extinc- 
tion in the NGC 3310 starburst regions to be in the range 
0.18 < E{B - V^) < 0.23, based on the slope of the UV 
SED and published Balmer decrement measurements (P93) , 
respectively, and P93 estimated 0.13 < E{B - V) < 0.37 
based on spectroscopic measurements (see also S96, DOO). 
Since these extinction estimates cover both the YSCs and 
the ISM, they are consistent with being upper limits to the 
extinction towards individual YSCs. 

For their large-scale complexes, E02 assumed a mean 
age of ~ lO'^ yr, which is somewhat underestimated com- 
pared to our mean age of the individual clusters in these 
complexes of ~ 4 x 10^ yr. For the "Jumbo" region (E02 
complex 108-1-109), P93 derived an age of 1.45 x 10^ yr from 
spectrophotometry, which is consistent with the median age 



we derive for the clusters in this region, fjumbo ~ 1-3 x 10^ 

yr- 

Although the match between our individual cluster age 
estimates and the E02 results is good for clusters 8 (E02 
estimate ~ 20 — 30 Myr; our estimate ~ 26 Myr), 19 (~ 
20 - 30 vs. ~ 21 Myr) and 64 (10 - 50 vs. 33 Myr), and 
reasonable for clusters 36 (< 5 — 10 vs. ~ 9 — 10 Myr) and 
90 (10 - 50 vs. 14 Myr), it is poor for clusters 37 (~ 1 - 5 
vs. ~ 16 Myr; this is no surprise, because our models do not 
extend to such young ages), 49 (10 — 50 vs. 94 Myr), and 59 
(10 — 50 vs. 320 Myr). Owing to the strong age-dependence 
of the M/L ratio at those young ages, small differences in 
our age estimates will have significant consequences for the 
final mass determinations, for ages in the range determined 
for the NGC 3310 clusters. In addition, the difference in the 
assumed distance modulus to NGC 3310 between ourselves 
{D = 13 Mpc; Sect. 2.1) and E02 {D = 18.7 Mpc) of ~ -0.8 
mag introduces an additional uncertainty to the clusters' 
luminosity of AlogL ~ 0.3, corresponding to a factor of 
~ 2 in luminosity, and therefore in the mass (via the M/L 
ratio) as well, in the sense that we expect the E02 mass 
estimates to be ~ 2x higher. 

Thus, from a comparison between our final masses and 
those of E02, we conclude that the factor of 2 — 3 between 
most of our mass estimates, in the sense that E02's masses 
are the higher, is not too surprising in view of the many un- 
certainties involved in (i) the difference in distance to NGC 
3310 assumed, (ii) the different methods used for our mass 
estimates, (iii) the different ways in which we obtained our 
aperture photometry and the subsequent photometric cali- 
brations (E02 converted their measurements to the standard 
ground-based system, while we continued in the STMAG 
system), and (iv) the different models employed (cf. Fritze- 
V. Alvensleben 2000) . Regarding the latter, we point out that 
these differences may have been augmented by our inclusion 
of nebular emission (which we have shown to be important 
for ages younger than ~ 3 x 10^ yr; Anders et al. 2002, 
Anders & Fritze-v. Alvensleben 2003) and our approach to 
solve for all of our free parameters simultaneously. 

This comparison clearly shows that while we can obtain 
robust relative distributions of masses, ages, and metallic- 
ities for a given star cluster system fairly easily (as shown 
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Table 3. Comparison of our derived cluster parameters with those of the E02 SSC candidates. 
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Figure 6. (top panel) Age dependence of the model SEDs for 
solar metallicity. (bottom panel) Metallicity dependence of the 8 
Myr-old model SED. we also show the effects of adding E(i3 — 
V) = 0.5 mag extinction to the solar-metallicity SED (dashed 
line with solid bullets). The HST filters are spaced according to 
their central wavelengths. The grey vertical lines are drawn at the 
central wavelengths of the F606W (or F555W) and F814W HST 
filters often used for age dating of individual clusters. 
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Figure 7. Age distribution of the NGC 3310 star clusters based 
on a comparison of their SEDs with the Starburst99 models and 
fixed, solar metallicity. The dashed Gaussian distribution is our 
best-fit age distribution from a comparison with the Gottingen 
SSP models, leaving the metallicity as a free parameter. The error 
bars represent Poissonian uncertainties in the numbers of clusters 
in each age bin. 



throughout this paper), the absolute cahbration, and there- 
fore the absolute mass, age and metallicity scales can only be 
fixed robustly if we have access to independent spectroscopic 
measurements of these parameters. 



6.3 The NGC 3310 Cluster Luminosity Function 

The importance of correcting the observed CLFs to a com- 
mon age cannot be overemphasized. Age spread effects in 
young cluster systems, and therefore the combined effects 
of ongoing cluster formation, evolutionary fading, and the 
onset of cluster disruption, affect the observed CLF (cf. 
Meurer 1995, Fritze-v. Alvensleben 1999, de Grijs et al. 
2001, 2003a,b). This implies that the CLF observed in such 
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a system represents merely a temporal situation, rather than 
a characteristic property of a coeval cluster system. 

Using the age estimates obtained for the individual clus- 
ters, we corrected the F606W magnitudes of our cluster sam- 
ple to a common age of lO'^"'' yr using the Anders & Fritze- 
V. Alvensleben (2003) models properly folded through the 
HST/WFPG2 F606W filter response curve, for the "best-fit" 
age distribution of Fig. 5a. In order to provide a comparison 
with the CLFs of other young cluster systems, wo employed 
the Lcvenberg-Marquardt method (Marquardt 1963) to fit 
power-law slopes of the form N{L) oc L" (where L is the 
luminosity of a cluster, and a the power-law slope) to both 
CLFs. 

In the F-band luminosity range 10*' < Lp606w/i© < 
lO'^ (17.7 < F606W < 20.2 mag), the slope of the CLF is 
a = —1.4 ±0.2, and a = —1.8 ±0.4 for the full, uncorrected 
sample of clusters, and for the "best-fit" sample corrected 
to a common age of log(age/yr) — 7.5, respectively. 

In de Grijs et al. (2001) and Parmentier et al. (2002) 
we showed that the constant-age F-band CLF for the bright 
clusters in M82's post-starburst region "B" roughly follows a 
power law with a slope in the range —1.4 < a < —1-2, which 
is consistent with the power-law slopes of other young CLFs, 
and within the uncertainties also with the NGC 3310 CLF. 
On the other hand, DOO found a significantly shallower slope 
of a ~ —1 for the Ha-bright circumnucleax star forming 
regions in NGC 3310 (but we note that these are likely the 
youngest objects populating the CLF, with ages < 6 — 12 
Myr, depending on metallicity), while E02 obtained slopes 
as steep as (-2.4 ± 0.04) < a < (-2.2 ± 0.03) from their 
B and if-band cluster measurements, which is - within the 
uncertainties - consistent with our V-band determination. 
The main uncertainties that arc expected to affect the CLF 
slope in this latter case are the luminosity range over which 
the slopes were obtained, and the fact that E02 did not 
correct their CLFs to a common age before measuring their 
slopes. 



timisation technique to simultaneously determine the best 
combination of age, extinction and metallicity. 

Our main results can be summarised as follows. For a 
star cluster system as young as that in NGC 3310, 

• the peak of the age distribution is robustly reproduced 
for all of our choices of passband combinations, except for 
the red-only selected filter combination; 

• red-dominated passband combinations result in signif- 
icantly different (uncertain) age solutions, due to the weak 
time dependence of the NIR magnitudes and due to unavoid- 
able ambiguities in the modelling of the thermally pulsing 
AGB phase; they produce a significant wing of older clusters. 

• blue-selected passband combinations tend to result in 
age estimates that are slightly skewed towards younger ages, 
compared to passband combinations that also include red- 
der passbands. This is due to the combination of observar 
tional selection effects and to the age-metallicity degener- 
acy, which is also of importance in passband combinations 
containing only optical filters. 

• for ages 6 < log(age/yr) < 9, we conclude that if one 
can only obtain partial coverage of a star cluster's SED, 
an optical passband combination of at least 4 filters includ- 
ing both blue and red optical passbands results in the best 
balanced and most representative age distribution. While 
blue-selected passband combinations lead to age distribu- 
tions that are slightly biased towards younger ages, red and 
in particular NIR-dominated passband combinations should 
be avoided. 

• The physical effect most limiting the accuracy of our 
age determinations, for ages up to ~ 10*^ yr, is the age- 
metallicity degeneracy. This is clearly illustrated in our com- 
parison with results obtained assuming a fixed, solar metal- 
licity. If the actual cluster metallicities are significantly sub- 
solar, this assumption results in a markedly different age 
distribution, biEised towards younger ages, compared to the 
distribution obtained by leaving the metallicity as a free pa- 
rameter. 



7 SUMMARY AND CONCLUSIONS 
7.1 Systematic uncertainties 

In this paper we have investigated the systematic uncer- 
tainties involved in using broad-band UV, optical and NIR 
observations to derive age, mass, metallicity and extinction 
distributions for extragalactic star cluster systems. The large 
majority of extragalactic star cluster studies done to date 
have essentially used two or thrcc-iiassbaud aperture pho- 
tometry, combined with theoretical stellar population syn- 
thesis models to obtain age estimates. The accuracy to which 
this can be done obviously depends on the number of differ- 
ent (broad-band) filters available as well as, crucially, on the 
actual wavelengths and the wavelength range covered by the 
observations. Understanding these systematic uncertainties 
is therefore of the utmost importance for a well-balanced 
interpretation of the properties of such star cluster systems. 

We focused our analysis on the nearby, well-studied 
starburst galaxy NGC 3310, known to harbour large num- 
bers of young star clusters, for which we obtained multi- 
passband observations from the UV to the NIR from the 
HST Data Archive. We used our evolutionary synthesis op- 



7.2 The NGC 3310 cluster system 

NGC 3310 is a local, very active starburst gala^xy with high 
global star formation efficiency, which has experienced very 
recent star formation in star clusters and Hll regions, in the 
last < 10'^ - 10* yr. In the inner FoV covered by all of our 
passbands, we detected some 300 star cluster candidates, 
while an additional ~ 100 clusters were detected in the outer 
FoV covered by the WFPC2 filters only, all of them affected 
by only moderate levels of NGC 3310- internal extinction to- 
wards the clusters, E{B—V) < 0.1 mag. The age distribution 
derived for these clusters indicates that NGC 3310 under- 
went a significant burst of cluster formation some 3 x 10^ yr 
ago. The actual duration of the burst of cluster formation, 
which we estimate to have lasted for almost 20 Myr, may 
have been shorter because uncertainties in the age determi- 
nations may have broadened the peak. It appears, therefore, 
that the peak of cluster formation in NGC 3310 coincides 
closely with the suspected galactic cannibalism or last tidal 
interaction. 

The clusters older than ~ lO'^ yr are smoothly dis- 
tributed throughout the galactic centre. However, the 
younger clusters, with ages log(age/yr) < 7.1, are predomi- 
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nantly concentrated in the Jumbo region and in the northern 
spiral arm. This is consistent with the very young ages in- 
dependently derived for these regions, t < 10 Myr, although 
both regions also contain clusters spanning the entire age 
range observed for the NGC 3310 star cluster system. 

The CLF is the result of a complex interplay of the 
intrinsic cluster mass distribution, age spread and cluster 
disruption processes. We therefore determined the individ- 
ual cluster masses by scaling our model SEDs for the best- 
fit age, extinction values and mctallicity estimates to the 
observed SEDs. We estimate the cluster system to have a 
median mass of (log(m/M0)) ~ 5.25 ± 0.1, not including 
systematic uncertainties introduced by the uncertainties in 
the low-mass IMF slope. 

Our metallicity determinations are strongly dominated 
by (significantly) subsolar metallicities, which is consistent 
with independent metallicity measurements. There is some 
evidence that the most actively star forming regions, in par- 
ticular the Jumbo region and the northern spiral arm, arc 
predominantly composed of lower-abundance star clusters. 
The V-band CLF slope in the range 10^ < Lpeoew/.i'© < 
10^ for the NGC 3310 star cluster system is apeoew ~ 
— 1.8 ± 0.4, which is consistent with the power-law slopes 
of other young CLFs. 

Finally, we point out that our estimates for the ages, 
masses, metallicities and extinction values of the YSCs in 
NGC 3310 closely match previous, independent determina- 
tions of those parameters, where available (e.g., Grothous 
& Schmidt-Kaler 1991, P93, Meurer et al. 1995, DOO, E02). 
This implies, to first order, that this type of analysis is rea- 
sonably robust. 



ACKNOWLEDGMENTS 

This paper is based on new and archival observations with 
the NASA/ESA Hubble Space Telescope, obtained at the 
Space Telescope Science Institute (STScI), which is operated 
by the Association of Universities for Research in Astronomy 
(AURA), Inc., under NASA contract NAS 5-26555. This 
paper is also partially based on ASTROVIRTEL research 
support, a project funded by the European Commission un- 
der 5FP Contract HPRI-CT-1999-00081. This research has 
made use of NASA's Astrophysics Data System Abstract 
Service. PA is partially supported by DFG grant Fr 916/11- 
1; PA also acknowledges partial funding from the Marie 
Curio Fellowship programme EARASTARGAL "The Evolu- 
tion of Stars and Galax;ies" , funded by the European Com- 
mission under 5FP contract HPMT-CT-2000-00132. VAT 
and RAW acknowledge support from NASA grants GO- 
8645.*, awarded by STScI. RdG acknowledges preliminary 
analysis of the Cycle 9 HST UV data by Xsitaaz Chadcc as 
part of the 2001 PPARC/Cambridge International Under- 
graduate Summer School. 



REFERENCES 

Anders P., Fritze-v. Alvensleben U., 2003, A&A, in press (astro- 
ph/0302146) 

Anders P., Pritze-v. Alvensleben U., de Grijs R., 2002, in: The 



Evolution of Galaxies. Ill - From Simple Approaches to Self- 
Consistent Models, cds. Hensler G., et al., Kluwer: Dordrecht, 
in press (astro-ph/0210223) 
Anders P., Pritze-v. Alvensleben U., de Grijs R., 2003, MNRAS, 
in prep. 

Athanassoula E., 1992, MNRAS, 259, 345 

Balick B., Heckman T., 1981, A&A, 96, 271 (BH81) 

Bertola P., Sharp N.A., 1984, MNRAS, 207, 47 

Bik A., Lamcrs H.J.G.L.M., Bastian N., Panagia N., Romaniello 

M., 2003, A&A, 397, 473 
Biretta J. A., ct al., 2000, WFPC2 Instrument Handbook, Version 

5.0 (Baltimore; STScI) 
Bottinelli L., Gougenheim L., Paturel G., de Vaucouleurs G., 

1984, A&AS, 56, 381 
Boutloukos S.G., Lamers H.J.G.L.M., 2003, MNRAS, 338, 717 
Braine J., Combes P., Casoli P., Dupraz C, Gerin M., Klein U., 

Wielebinski R., Brouillet N., 1993, A&AS, 97, 887 
Bruzual G., Chariot S., 2000, updated version of Bruzual G., 

Chariot S., 1996, in: Lcitliercr, C, et al. 1996, PASP, 108, 

996 (AAS CDROM Series 7) (BCOO) 
Calzetti D., Armus L., Bohlin R.C., Kinney A.L., Koornneef J., 

Storchi-Bergmann T., 2000, ApJ, 533, 682 
Calzetti D., Kinney A.L., Storchi-Bergmann T., 1994, ApJ, 429, 

582 

Carlson M.N., et al., 1998, AJ, 115, 1778 
Code A., Welch G., 1982, ApJ, 256, 1 

Conselice C.J., Gallagher J.S., Calzetti D., Homeier N., Kinney 

A.L., 2000, AJ, 119, 79 
de Grijs R., Bastian N., Lamers H.J.G.L.M., 2003a, MNRAS, in 

press (astro-ph/0211420) 
de Grijs R., Bastian N., Lamers H.J.G.L.M., 2003b, ApJ, 583, 

L17 

de Grijs R., Johnson R.A., Gilmore G.P., Prayn CM., 2002, MN- 
RAS, 331, 228 

de Grijs R., O'Connell R.W., Gallagher J.S., 2001, AJ, 121, 768 
Diaz A. I., Alvarez Alvarez M., Tcrlevicli E., Tcrlevicli R., Sanchez 

Portal M., Aretxaga I., 2000, MNRAS, 311, 120 (DOO) 
Elmegrecn B.C., Efremov Y.N., 1997, ApJ, 480, 235 
Elmegreen D.M., Chromey P.R., McGrath E.J., Ostenson J.M., 

2002, AJ, 123, 1381 (E02) 
Elson R.A.W., Pall S.M., 1985, PASP, 97, 692 
Eskridge P.B., et al., 2003, ApJ, in press (astro-ph/0211494) 
Pabbiano G., Peigelson E., Zamorani G., 1982, ApJ, 256, 397 
Pall S.M., Rees M.J., 1977, MNRAS, 181, 37P 
Pall S.M., Rces M.J., 1985, ApJ, 298, 18 
Pall S.M., Zhang Q., 2001, ApJ, 561, 751 

Ferguson A.M.N., Gallagher J.S., Wyse R.P.G., 1998, AJ, 116, 
673 

Pritze-v. Alvensleben U., 1998, A&A, 336, 83 
Pritze-v. Alvensleben U., 1999, A&A, 342, L25 
Pritze-v. Alvensleben U., 2000, in Langon A., Boily C, eds., ASP 
Conf. Ser. Vol. 211, Massive Stellar Clusters, Astron. Soc. 

Pac, San Francisco, p. 3 
Fritzc-v. Alvensleben U., Gerhard O.E., 1994, A&A, 285, 751 
Gnedin O.Y., Ostrikcr J. P., 1997, ApJ, 474, 223 
Grothcus H.G., Schmidt-Kaler T., 1991, A&A, 242, 357 
Harris W.E., 1996, AJ, 112, 1487 

Harris W.E., 2001, in: Star Clusters, Saas-Pee Advanced Course 

28, Springer- Verlag, 223 
Harris W.E., Harris G.L.H., McLaughlin D.E., 1998, AJ, 115, 

1801 

Harris W.E., Pudritz R.E., 1994, ApJ, 429, 177 

Heckman T.M., Balick B., 1980, A&AS, 40, 295 

Hibbard J.E., Milios J.C., 1995, AJ, 110, 140 

Kinney A.L., Bohlin R.C., Calzetti D., Panagia N., Wyse R.F.G., 

1993, ApJS, 86, 5 
Kobulnicky H.A., Skillman E.D., 1995, ApJ, 454, L121 
Kregel M., Sancisi R., 2001, A&A, 376, 59 (KSOl) 



© 2002 RAS, MNRAS 000, 000-000 



Star Clusters in Starburst Galaxies. I. 17 



Kroupa P., Tout C.A., Gilmorc G.F., 1993, MNRAS, 262, 545 
Lcithcrcr C, Scliacrcr D., Goldadcr J.D., Gonzalez Delgado R.M., 

Robert C, Kune D.F., de Mello D.F., Devest D., Heckman 

T.M., 1999, ApJS, 123, 3 (Starburst99) 
Maoz D., Bajrth A. J., Ho L.C., Sternberg A., Filippenko A.V., 

2001, AJ, 121, 3048 
Maoz D., Barth A. J., Sternberg A., Filippenko A.V., Ho L.C., 

Macchetto F.D., Rix H.-W., Schneider D.P., 1996, AJ, 111, 

2248 

Marquardt D.W., 1963, J. of the Soc. for Industrial and Applied 

Mathematics, 11, 431 
Martin C.L., Kennicutt R.C., Jr., 1995, ApJ, 447, 171 
McLaughlin D.E., Pudritz R.E., 1996, ApJ, 457, 578 
Mengel S., Lehnert M.D., Thatte N., Tacconi-Garman L.E., Gen- 

zel R., 2001, ApJ, 550, 280 
Meurer G.R., 1995, Nat, 375, 742 

Meurer G.R., Heckman T.M., Leitherer C., Kinney A., Robert 

C., Garnctt D.R., 1995, AJ, 110, 2665 
Miller B.W., Whitmorc B.C., Schwcizcr F., Fall S.M., 1997, AJ, 

114, 2381 

Mirabel I.F., ct al., 1998, A&A, .333, LI 
Murali C., Weinberg M.D., 1997a, MNRAS, 288, 749 
Murah C., Weinberg M.D., 1997b, MNRAS, 288, 767 
Murali C., Weinberg M.D., 1997c, MNRAS, 291, 717 
Mulder P.S., van Driel W., 1996, A&A, 309, 403 {MvD96) 
Mulder P.S., van Driel W., Braine J., 1995, A&A, 300, 687 
(MvDB95) 

Oey M.S., Kennicutt R.C., Jr., 1993, ApJ, 411, 137 

Ostrikcr J. P., Gncdin O.Y., 1997, ApJ, 487, 667 

Parmenticr G., dc Grijs R., Gilmore G., 2002, MNRAS, submitted 

Pastoriza M.G., Dottori H.A., Tcrlcvich E., Terlevich R., Diaz 

A.I., 1993, MNRAS, 260, 177 (P93) 
Finer B.G., Stone J.M., Teuben P.J., 1995, ApJ, 449, 508 
Puxley P.J., Hawarden T.G., Mountain C.M., 1990, ApJ, 364, 77 
Richer M.G., McCall M.L., 1995, ApJ, 445, 642 
Ratnatunga K.U., Bahcall J.N., 1985, ApJS, 59, 63 



Schlcgcl D.J., Finkbeiner D.P., Davis M., 1998, ApJ, 500, 525 
Schulz J., Fritze-v. Alvensleben U., MoUer C.S., Fricke K.J., 2002, 
A&A, 392, 1 

Schweizer F., Miller B.W., Whitmore B.C., Fall, S.M., 1996, AJ, 
112, 1839 

Schweizer F., Seitzer P., 1988, ApJ, 328, 88 

Scuderi S., Panagia N., Gilmozzi R., Challis P.M., Kirshner R.P., 

1996, ApJ, 465, 956 
Smith D.A., ct al., 1996, ApJ, 473, L21 (S96) 
Stetson P.B., 1987, PASP, 99, 91 
Telesco CM., Gatley L, 1984, ApJ, 284, 557 (TG84) 
Terlevich E., Diaz A.L, Pastoriza M.G., Terlevich R., Dottori H., 

1990, MNRAS, 242, 48P 
van der Kruit P.C., 1976, A&A, 49, 161 
van der Kruit P.C., de Bruyn A.G., 1976, A&A, 48, 373 
van Zee L., Salzcr J. J., Haynes M.P., O'Donoghue A. A., Balonek 

T.J., 1998, AJ, 116, 2805 
Vcspcrini E., 2000, MNRAS, 318, 841 
Vesperini E., 2001, MNRAS, 322, 247 
Walker M.F., Chincarini G., 1967, ApJ, 147, 416 
Whitmore B.C., Schweizer F., 1995, AJ, 109, 960 
Whitmore B.C., Schweizer F., Leitherer C, Borne K., Robert C, 

1993, AJ, 106, 1354 
Whitmore B.C., Schweizer F., Kundu A., Miller B.W., 2002, AJ, 

124, 147 

Whitmore B.C., Sparks W.B., Lucas R.A., Macchetto F.D., 

Birctta J.A., 1995, ApJ, 454, L73 
Windhorst R.A., et al., 2001, NASA press release, STScI-PROl- 

04, 11 January 2001 

(http: / /hubblcsitc.org/ncwsccntcr/archivc/2001/04/) 
Windhorst R.A., ct al., 2002, ApJS, 143, 113 
Zaritsky D., Kennicutt R.C., Jr., Huchra J.P., 1994, ApJ, 420, 87 
Zepf S.E., Ashman K.M., English J., Freeman K.C., Sharpies 

R.M., 1999, AJ, 118, 752 
Zezas A.L., Georgantopoulos I, Ward M.J., 1998, MNRAS, 301, 

915 



© 2002 RAS, MNRAS 000, 000-000 



